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ABSTRACT 

The  objectives  of  the  present  report  are  to  understand  thermal  injury 
and  to  provide  methods  for  evaluation  of  skin  damage  due  to  thermal 
exposures.     This  report  reviews  the  literature  of  thermal  injury  for  the 
period  from  the  late  1940 's  to  the  1970' s.     Examination  of  the  literature 
into  the  1980 's  found  that  the  quantification  of  thermal  injury  is  still 
based  on  the  earlier  work.     This  report  covers  a  mathematical  and 
experimental  analysis  of  heat  transfer  to  the  skin,   the  studies  of  thermal 
properties  of  the  skin  and  its  physiological  reactions  to  burns.  Methods 
to  estimate  skin  damage  incurred  during  a  burn  were  also  included. 
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1 .  INTRODUCTION 

Thermal  injury  was  a  major  source  of  morbidity  and  mortality  during 
World  War  II.     Since  that  time,  attentions  have  been  paid  to  the  studies  of 
skin  damage  due  to  thermal  exposure  in  order  to  develop  protective  devices 
and  to  give  medical  care  to  the  military.     Consequently,  numerous  reports 
on  this  topic  appeared.     Today,  statistics  indicate  that  thermal  injury  is 
one  of  the  major  types  of  trauma  in  the  modern  society.     Thermal  injury 
frequently  occurs  at  industrial  sites  where  the  risk  of  fires  is  high. 
Burns  from  clothing  fires  also  occur  most  often  in  the  home. 

One  method  to  minimize  bum  injury  is  to  use  heat  resistant  fabrics 
for  safety  apparel.     Instrumented  mannikins  are  usually  employed  in 
experiments  to  evaluate  protective  clothing.     From  the  measurements  of  heat 
transfer  through  fabrics  to  the  mannikln,  the  severity  of  burns  can  be 
assessed  by  means  of  the  available  knowledge  regarding  the  relationship 
between  the  heat  transfer  and  Its  resultant  skin  damage. 

Thermal  damage  to  the  skin  occurs  as  a  result  of  an  increase  in  tissue 
temperature  above  an  injurious  level  within  a  finite  period  of  time.  The 
severity  of  bums  depends  on  the  intensity  and  duration  of  exposures.  The 
mode  of  heat  transfer  can  be  conduction,  convection,  radiation,  or  any 
combination  of  the  above  mechanisms. 

To  understand  the  problem  of  thermal  injury  and  to  provide  adequate 
information  for  assessment  of  skin  damage  due  to  thermal  exposure  through 
the  use  of  an  instrumented  mannikin,  the  present  report  reviews  the 
literature  of  thermal  injury.    The  pioneering  works  of  many  investigators 
(around  1950  to  1970)  will  be  presented.    A  mathematical  and  physical 
analysis  of  the  physiological  response  of  skin  to  burns  will  be  covered. 
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Chapter  2  briefly  describes  the  anatomy  of  human  skin.     Chapter  3 
illustrates  the  response  of  skin  temperature  to  various  conditions  of  heat 
exposure.     Mathematical  expressions  describing  the  change  of  skin  tempera- 
ture under  different  conditions  are  deduced.     The  thermal  properties  of  the 
skin  are  discussed  in  Chapter  4.     The  variations  of  these  properties  with 
layer,  site,  moisture  content,  irradiation  and  temperature  are  explained. 
In  Chapter  5,   the  quantitative,  qualitative  and  sensational  effects  of  heat 
on  the  skin  are  studied.     Information  concerning  the  thresholds  of  burns, 
blistering  and  pain  as  related  to  the  intensity  and  durations  of  thermal 
exposure  is  given.     Factors  that  may  alter  the  vulnerability  of  skin  to 
thermal  injury  are  also  discussed.     Chapter  6  summarizes  the  physical  and 
chemical  processes  that  participate  in  cell  death  during  a  burn.     A  rate 
process  model  is  developed  in  this  chapter  to  predict  the  skin  damage. 
Relationship  between  pain  intensity  and  skin  damage  is  shown.     Chapter  7 
demonstrates  some  of  the  numerical  procedures  that  are  available  for 
determination  of  the  severity  of  bums  from  experimental  data  using  digital 
computers .     This  is  followed  by  a  brief  sunmiary  of  the  review  in  Chapter  8 . 
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2.       THE  HUMAN  SKIN 

The  skin  covers  the  entire  body  and  is,  therefore,  as  large  as  the 
body  itself.     For  most  adults,  the  area  is  about  1.6  to  1.9  square  meters 
[1] .     Table  1  shows  the  weighting  of  area  for  each  portion  of  the  body  [2] . 
The  functions  of  the  skin  are  diverse  and  include  protection,  excretion, 
sensation,  maintaining  fluid  and  electrolyte  balance  and  body  temperature. 

TABLE  1 

Weighting  of  Areas  with  Respect  to  Total  Skin  Surface, 
From  Stoll  and  Chianta  [2] 

%  of  Total  Area 


Head  7 

Trunk  35 

Arms  14 

Hands  5 

Thighs  19 

Legs  13 

Feet  7 


2.1    The  Structure  of  Skin 

The  skin  is  made  up  of  two  main  layers:     an  outer  layer,  the 
epidermis,  and  an  inner  layer,  the  dermis  or  corium.     The  tissues  under  the 
skin,  such  as  the  fatty  tissue  and  muscle,  are  called  the  subcutaneous  or 
hypodermic  tissues.     The  tissues  within  the  layers  of  the  skin  are  called 
the  intracutaneous  or  cutaneous  tissues.    The  structures  of  the  epidermis 
and  the  dermis  are  depicted  in  Figures  1  and  2,  respectively. 

The  epidermis  consists  of  squamous  epithelial  tissue  and  is  divided 
into  characteristic  strata  (Figure  1) .    The  stratum  germinativum  or  the 
basal  layer  is  the  deepest  stratum  from  which  new  cells  are  produced.  The 
cells  move  upward  as  they  form  and  gradually  flatten.     In  the  stratum 
granulosum,  the  cells  die  and  their  nuclei  degenerate.    The  dead  cells  are 
converted  into  a  tough  fibrous  protein  called  keratin  which  is  water- 


Figure  2      The  Dermis,  from  Gardner  and  Osburn  [3] 


repellant  and  resistant  to  abrasion.  The  keratin  is  formed  in  the  stratum 
corneum  (the  horny  layer)  in  which  the  new  keratin  keeps  pushing  upward  to 
replace  the  old  ones  that  continually  flake  off  from  the  skin  surface. 

The  dermis  is  composed  of  fibrous  connective  tissue  which  provides 
strength  and  elasticity  for  the  skin  (Figure  2).     The  glands  of  the  skin, 
the  hair  follicles,  and  the  blood  vessels  are  all  located  within  this 
layer.     The  dermis  is  further  divided  into  a  superficial  papillary  layer 
and  a  deeper  reticular  layer.     The  papillary  layer  presses  against  the 
epidermis  so  that  the  blood  vessels  in  the  papillary  layer  are  in  close 
contact  with  the  epidermis.     As  a  result,  nutriments  and  oxygen  reach  the 
epidermis  by  diffusion.     The  papillary  layer  is  also  thrown  up  into 
papillary  ridges,  especially  in  the  fingers,  palms  and  soles.  These  ridges, 
covered  with  epidermis ,  form  an  unique  pattern  on  the  skin  surface  in  each 
individual.     These  patterns  of  the  ridges  are  known  as  finger  prints,  palm 
prints  or  sole  prints. 

The  root  of  the  hair  is  embedded  in  the  reticular  layer.  Together 
with  its  coverings  (a  continuation  of  the  stratum  germinativum) ,  the  root 
forms  the  hair  follicle.    At  the  bottom  of  the  follicle  is  a  loop  of  capil- 
laries enclosed  in  a  connective  tissue  papilla.    The  epithelial  cells  over 
the  papilla  reproduce  and  undergo  keratinization  (a  process  of  keratin 
formation)  to  form  the  hair  shaft.    A  small  band  of  involuntary  muscle 
fibers,  the  arrector  pili,  connects  the  hair  follicle  to  the  papillary 
layer.    When  these  muscles  contract  in  response  to  cold  or  fright,  they 
pull  the  hair  erect.    As  a  result,  the  skin  around  the  hair  is  raised  and 
goose  pimples  are  formed.    Two  or  more  sebaceous  glands  cluster  beside  each 
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hair  follicle.     These  glands  secrete  an  oily  product,  sebum,  to  keep  the 
hairs  soft  and  pliable. 

Sweat  glands  are  long  tubular  glands.     The  coiled  secretory  portions 
of  the  glands  are  located  in  the  reticular  layer  or  even  in  the  sub- 
cutaneous layer.     Sweat  secretion  eliminates  water  from  the  body  to 
maintain  fluid  and  electrolyte  balance  and  body  temperature. 

Millions  of  sensory  receptors  are  distributed  throughout  the  skin. 
These  receptors  are  the  distal  ends  of  dendrites  of  all  sensory  nerve  cells 
(neurons) .     The  simplest  receptors  are  merely  the  free  endings  of  sensory 
dendrites.     They  are  mainly  in  the  dermis.     Some  of  them  surround  the  root 
of  the  hair  follicle  in  linear  and  circular  fashion,  and  some  extend  to  the 
epidermis.     The  free  nerve  endings  are  those  which  mediate  touch,  coldness, 
hotness  and  pain  such  as  that  caused  by  thermal  injury.     A  more  complete 
description  of  the  structure  of  the  skin  and  its  sensory  organ  is  in 
References  [1]  and  [3]. 

2,2    The  Thickness  of  Skin 

Unlike  the  skin  structure,  the  thickness  of  the  skin  varies 
considerably  with  age,  sex,  race  and  site.     Southwood  [4]  studied  the 
variations  of  the  thickness  of  both  epidermis  and  dermis  in  various  regions 
from  both  sexes  and  at  different  ages.     He  commented  that  the  thickness  of 
the  epidermis  was  inversely  proportional  to  the  hair  density.    Thus  the 
areas  with  the  thickest  epidermis  were  those  having  no  hairs,  for  example: 
the  palms  and  the  soles.     In  his  paper,  the  only  comment  on  the  difference 
in  skin  thickness  between  races  is  that  the  epidermis  of  Negro  was  thicker 
than  that  of  the  European.    His  results  are  shown  in  Tables  2  and  3. 


Table  2      The  Thickness  of  the  Epidermis  in  fxm,   from  Southwood  [4] 
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Table  3      The  Thickness  of  the  Dermis  in  /xm,  from  Southwood  [4] 
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As  shown,  men  have  thicker  dermis  than  women.     The  epidermis  is 
thickest  in  the  sole.     But  the  dermis  is  thickest  in  the  trunk.     At  birth 
and  during  childhood,  the  epidermis  is  relatively  thin.     It  becomes  thicker 
at  about  puberty  and  maintains  this  thickness  until  50  to  60  years  old. 
Then  it  becomes  thinner  again.     However,  the  thickness  of  the  dermis 
increases  with  age  to  a  maximum  at  about  30  to  40  years  old  and  sub- 
sequently decreases.     In  general,  the  thickness  of  the  epidermis  and  dermis 
ranges  from  20  to  1400  /xm  and  400  to  2500  ^m,  respectively. 


3.       CHANGE  OF  SKIN  TEMPERATURE  IN  RESPONSE  TO  VARIOUS  CONDITIONS  OF  HEAT 
EXPOSURE 

When  a  body  is  exposed  to  a  heat  flux,   the  temperature  of  the  skin 
increases.     If  the  heat  flux  is  intense  enough  to  raise  the  skin  tempera- 
ture beyond  a  safety  limit  within  a  period  of  time,   thermal  injury 
(cutaneous  burn)  occurs.     To  prevent  the  occurrence  of  thermal  injury  by 
protective  devices,  a  knowledge  of  the  responses  of  skin  temperature  to 
different  kinds  of  heat  applications  is  essential.     In  the  past,  mathema- 
tical expressions  have  been  formulated  to  describe  the  relationships 
between  the  heat  transfers  and  the  skin  temperature  changes .  Experimental 
results  of  the  responses  of  skin  temperature  to  different  kinds  of  exposure 
have  also  been  determined. 

3 . 1    Mathematical  Deductions 

Although  there  are  numerous  complex  factors  that  are  involved  in  the 
transfer  of  heat  to  and  through  the  skin  such  as  variations  of  skin  thick- 
ness, existing  temperature  gradients  within  the  skin,  blood  flow  rate 
through  the  various  skin  layers,  density,  heat  capacity  and  therm il  con- 
ductivity of  the  different  layers,  Henriques  and  Moritz  [5]  mentioned  that  * 
it  was  possible  to  derive  a  first  approximation  of  the  time- temperature 
relationship  in  the  basal  epidermal  layer  or  within  the  epidermis. 
Buettner  [6]  summarized  the  mathematical  deductions  for  various  conditions 
of  heat  exposure,  including  a  hot  surrounding  medium,  direct  thermal 
contact,  non- penetrating  radiant  heat,  penetrating  heat,  and  others. 

They  considered  the  skin  as  a  solid  with  a  plane,  infinite  boundary 
and  assumed  its  material  properties  to  be  independent  of  depth.     Then  the 
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solutions  can  be  obtained  from  the  energy  balance  equation  of  a  skin 
element  which  is  simply  the  common  differential  equation  of  heat 
conduction: 

where  T  -  temperature,  "C 

t  -  time,  sec. 

X  —  distance  below  the  surface,  cm 
a  -  k/pc  -  diffusivity,  cm^/sec 
k  -  heat  conductivity,  cal/cm  sec  "C 
p  -  density,  g/cw? 

C  -  specific  heat  capacity,  cal/g  "C 

Equation  (1)  can  be  solved  by  many  methods  such  as  Fourier  Integral, 
Laplace  Transformation  and  similarity  solution,  provided  that  the  boundary 
and  initial  conditions  are  known  [7,  8,  9]. 
3.1.1  Hot  Surrounding  Medium 

In  this  case,  the  body  is  surrounded  completely  by  an  envelope  of 
ambient  hot  fluid  (e.g.  warm  air).    Assuming  a  zero  temperature  gradient 
within  the  skin  at  t  -  0,  the  initial  and  boundary  conditions  are: 
t  -  0,  X  >  0,  T  -  To  (2) 

t>0.  x-0,  -  MTa  -  Ts)  (3) 

Then  the  temperature  T  at  depth  x  and  time  t  by  solving  Eqs.   (1),   (2)  and 
(3)  is: 

T  -T  ,    ,.2^  ^  2 


a  -  ,    x    .    .     (q  h  t/k  +hx/k)   ,     _     ,  x        .  h  vat. ,  . 


where  Tq  -  initial  skin  temperature,  'C 
Ta  -  ambient  temperature,  "C 
Tg  -  skin  surface  temperature,  "C 

h    -  heat  transfer  coefficient  between  the  skin  surface  and  the 
surroundings ,  cal/cm^ . sec . " C 


erf  (u)  -  -7-  e"^    dy  -  Error  Function 

Vtt  o 

erfc  (u)  -  1  -  erf(u)  -  complementary  Error  Function 
u  -  upper  limit  of  the  Integral  of  the  Error  Function 
y  -  dummy  variable 

The  results  of  the  Error  Function  for  different  values  of  u  are 
tabulated  in  p.  495  of  Reference  [8],     The  heat  transfer  coefficient,  h, 
depends  on  the  transfer  mechanisms.     Equations  that  describe  the  heat 
transfers  by  convecton  (q^)  and  radiation  (q^)  have  been  given  in  [5]. 
These  equations  are: 

qe  -  0.0026  (T^  -  Tg)^'^^  (5) 
qr  -  6.56  X  10'^^  (Tr  -  T^)  (6) 
Equation  (5)  is  the  convection  transfer  due  to  an  air  stream  with  a 
velocity  of  about  1.6  km/hr  to  a  cylindrical  object  of  about  30  cm  in 
diameter,  where  q^  is  in  cal/min.cm^,  T^  and  Tg  are  in  "C.     Equation  (6)  is 
the  radiation  transfer,   in  cal/min.cm^,  to  a  body  of  temperature  Tg ,   in  K, 
completely  enclosed  in  a  large  box  of  temperature  Tj-,  in  K. 

The  heat  transfer  coefficient,  h,  for  both  convection  and  radiation 
can  be  obtained  by  combining  the  contributions  of  convection  and  radiation 
to  heat  transfer  coefficient  as  computed  by  Eqs.   (5)  and  (6). 
3.1.2  Direct  Thermal  Cot>tact 

When  heat  transfer  occurs  between  the  skin  and  a  hot  body  by  contact, 
e.g.  boiling  water,  the  skin  surface  is  brought  immediately  to,  and 
maintained  at  a  constant  temperature  equal  to  the  hot  body  temperature.  In 
this  case,  the  heat  transfer  coefficient,  h,  is  almost  infinite,  and  the 
boundary  condition  is: 
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Solving  Eqs,   (1),   (2)  and  (7),  the  solution  for  this  case  is: 
■^s  ■  ^  X 

^  erf  (tt^)  (8) 


T     -  T  2  Jot:' 

s  o 

Equation  (8)  can  also  be  obtained  by  substituting  h  -  ®  and        -  Tg  in 
Eq.  (4). 

3.1.3  Nonpenetrating  Heat  Flow 

3,1.3.1    Constant  Heat  Flow 

In  this  case,  a  constant  heat  flux  (Q  in  cal/cm^sec)  is  absorbed  by 
the  upper  surface  of  the  skin,  for  example,  from  a  hot  radiator.  The 
boundary  condition  becomes: 

t  >  0,   X  -  0,    ai/ax  -  -  OA  (9) 

The  solution  of  Eqs .   (1),   (2)  and  (9)  is  then: 

T  -  T    -  2  ,  2  e-''^/^"'  -  X  [erfc  (:^)])  (10) 


o      k         "  IT  ^  ^2  Jot' 

3.1.3.2    Square -Wave  Pulse  of  Heat  Flow 

Equation  (10)  considers  only  the  heating  phase  during  thermal 
exposure.     If  the  decrease  of  temperature  after  the  heat  source  is  removed, 
are  taken  into  account.     Equation  (10)  can  be  extended  to  include  the 
cooling  phase.  For  a  square -wave  pulse  of  heat  source  described  as  follows: 
0<t<r,    Q-Qi"  constant 
t  <  0    and    t  >  r,     Q  -  0 
where  r  -  exposure  time,  sec. 

The  equation  suggested  by  Buettner  (given  in  Weaver  and  Stoll  [10])  is: 


^1  ,   ^    /a(t-r)     -xV4a(t-r)  ,     .  X 

-  r  ^  2  y  e  ^  -  X  (erfc       J  ^^^-r) 


)])  (11) 


14 

Equation  (11)  is  a  two-part  soluton.     During  the  heating  phase  0  <  t  <  r, 
the  first  term  on  the  right  hand  side  of  Eq.   (11)  is  real,  while  the  second 
term  is  imaginary  or  zero.     During  the  cooling  phase  t  >  r,  both  terms  are 
real.     As  time  approaches  infinity,  the  first  and  second  terms  approach 
zero  and  thus  T  approaches  Tq, 

3.1.4  Penetrating  Radiant  Heat 

Penetrating  radiant  heat  refers  to  the  radiation  from  the  sun  or  an 
atomic  bomb  blast.     Assuming  the  absorption  of  radiation  in  skin  layers  is 
analogous  to  many  familiar  types  of  radiation  absorption  such  as  that  in 
CO2  and  H2O,  the  decrease  in  intensity  of  the  penetrating  radiation 
transmitting  through  a  layer  is  proportional  to  the  thickness  of  the  layer 
and  the  intensity  of  the  radiation  at  that  point  [8].  Thus: 
dl  -  -  V  I  dx  (12) 
Where  I  —  penetrating  heat  flux,  cal/cm^sec 

1/  »  absorption  coefficient,  cm°^ 
Integrating  Eq.   (12)  from  the  surface  to  a  depth  x  gives: 

I  -  Is  e"*'''  (13) 
or 

dl/dx  -  -Igi/  e'*'^  (14) 
where  I^  -  heat  flux  impinging  at  the  surface.     In  this  case,  the  energy 
balance  equation  of  a  skin  element  is  no  longer  the  same  as  Eq.  (1).  The 
new  balance  equation  including  the  penetrating  effect  is: 

at  pC 


0 
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Nevertheless,   the  initial  condition  remains  unchanged  (Eq.    (2)),  but  the 
boundary  condition  becomes: 

OO.     x-0.  (16) 

where  R  -  non- penetrating  heat  flux  absorbed  at  the  surface  layer.  Then, 
the  solution  for  this  case  is 

I  e*'^"^ 

(17) 

I  e"*^^    2(1    +  R)    r-;r       2,.  x(I    +  R) 

s  ^  s  /  at     -X  /4at        ^  s        '  /     x  . 

-    +   \   V  —  e  -   5   erfc  (75 — 7==:) 

k»/  k  ir  k  2  Vat 

Numerical  results  from  the  above  equations  are  available  in  Refs.  5,  6 

and  10.     Tjrpical  time- temperature  relationshps  calculated  from  Equations 

(8),   (10),   (11)  and  (17)  are  given  in  Figures  3,  4,  5,  and  6  respectively. 

As  seen  in  Figure  6,  when  1/  is  greater  than  100  cm  ^,  the  increase  of 

surface  temperature  approaches  to  that  for  1/  -  «.     If  »/  -  «,  Equation  (17) 

reduces  to  Equation  (10),  the  non -penetrating  case. 


TCirCRATURE 


Figure  3      Temperature  Rises  at  Different 
Skin  Layers  Due  to  a  Direct 
Contact  Exposure.     Figures  on 
Curves  indicate  Time  in  Seconds. 
From  Buettner  [6] 


CONTACT  HEATING 
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Figure  4      Temperature  Rises  at  Different 
Skin  Layers  Due  to  a  Non- 
penetrating Constant  Heat 
Exposure.     Figures  on  Curves 
Indicate  Time  in  Seconds . 
From  Buettner  [6] 
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Figure  5      Temperature -Time  History  of 

Skin  Exposed  to  a  Square -Wave 
Pulse  of  Nonpenetrating 
Radiation.     From  Weaver  and 
Stoll  [10] 


Dote  3/23/66 

0  .  2  820 
A I  r  .0375 
A2  ..0360 

1  ..304 


'0    .25    50    .75   100  125  1.50   175  200 
T<ne  (Sec) 


17 


Figure  6      Increase  of  Skin  Surface 

Temperature  with  Time  During 
a  Penetrating  Radiation 
Exposure.     From  Buettner  [6] 
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3 . 2    Analysis  of  Theoretical  and  Experimental  Results 

Henriques  and  Moritz  [5]  computed  the  temperature  changes  at  the 
epidermal -dermal  interface  using  Equations  (4)  and  (8) .    They  assumed  that 
the  epidermal  thickness,  x,  was  80  /im  and  x/(2  Jol)  -  0.15  sec^'''^.  The 
computed  temperatures  for  the  skin  surface  exposed  to  a  hot  environment  and 
to  a  hot  body  by  direct  contact  are  shown  in  Tables  4  and  5  respectively. 

Comparing  the  data  in  the  tables,  one  finds  that  the  temperature  rise 
at  the  epidermal -dermal  junction  by  means  of  convection  and  radiation 
(Table  4)  is  much  slower  than  that  by  direct  contact  with  a  hot  body  (Table 
5).     In  the  latter  case,  the  temperatures  reach  steady  states  after  about 
30  to  60  sec  of  exposure.    Henriques  and  Moritz  concluded  that  "for  a  given 
source  temperature,  a  mechanism  that  enables  the  surface  temperature  to  be 
immediately  brought  to,  and  maintained  at,  the  source  temperature  has,  on  a 


time  basis,  at  least  a  thousand  times  greater  propensity  to  injure  the 
epidermis  than  a  heat  source  which  raises  the  skin  temperature  by  means  of 
radiation,  and  conduction  and  convection  of  relatively  immobile  air," 


Table  4      The  Computed  Time -Temperature  Relationships  at  the  Basal  Layer 
When  an  Entire  Animal  (~  30  cm  Diameter)  Is  Surrounded  by  an 
Envelope  of  Ambient  Air  and  Radiant  Heat  with  a  Constant 
Temperature  Source.     From  Henriques  and  Moritz  [5] 
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Table  5      The  Computed  Time -Temperature  Relationships  at  the  Basal  Layer 
When  the  Skin  Surface  Is  Immediately  Brought  to,  and  Maintained 
at,  a  Specific  Temperature.     From  Henriques  and  Moritz  [5] 
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Experiments  were  also  carried  out  by  Henriques  and  Moritz  to  study  the 
change  of  dermis -fat  interface  temperature  when  the  surface  of  porcine  skin 
was  immediately  brought  to,  and  maintained  at,  a  constant  temperature 
(Table  6).     After  10  min.  exposures  to  surface  temperatures  of  50  to  70°C, 
they  found  that  the  dermis  (corium)  thickness  increased  significantly  due 
to  the  excessive  accumulation  of  a  clear  yellowish  fluid  (edema)  in  the 
dermis.     When  the  surface  temperatures  were  maintained  below  50°C, 
accumulation  of  edema  fluid  did  not  appear.     Above  70 "C,  the  dermis  was 
denatured  so  that  these  mechanisms  were  destroyed.     The  presence  of  the 
edema  fluid  decreased  the  dermis  temperature  if  the  amount  of  edema  fluid 
presented,  more  than  compensated  for  the  energy  uptake.     This  effect  was 
maximal  when  the  skin  surface  was  maintained  at  about  60 "C. 


Table  6      The  Measured  Time -Temperature  Relationships  for  the  Dermal -Fat 
Interface  When  the  Skin  Surface  is  Immediately  Brought  to,  and 
Maintained  at,  a  Specific  Temperature.     From  Henriques  and 
Moritz  [5] 


Initial 

SUSM 

iff MC  tMipcfstarc  la  *C 

3$.o 

34.8 

34-8 

35» 

34-9 

34.3 

34.» 

34.5 

Daring 

65 

ezpotore 

45 

so 

5S 

60 

70 

80 

90 

X>«rak-Ut  laUHsec 

Uapcratm  la  *C 

34-8 

34-7 

34.$ 

34.6 

34.7 

34.« 

34^ 

M 

36 

3« 

39 

43 

46 

S« 

S3 

1^ 

o-S 

3S 

43 

45 

46 

5« 

6a 

65 

I.O 

S9-S 

4S 

47 

48 

S3 

65.6 

7« 

74 

».$ 

40 

47 

48 

47 

S3 

66.5 

7« 

77 

S.O 

40.S 

47 

49 

46 

54 

67 

74 

79 

4« 

47 

48.$ 

45 

^ 

67.S 

75 

79 

4* 

47 

47.S 

44.5 

58 

67.S 

77 

4S 

46.S 

47 

58 

XO.O 

4» 

47 

r 

49 

S9 

IhlrtaiM 

•feoriaab 

Initial 

s 

s 

s 

t 

t 

s 

At  tcnniaa- 

tionofez- 

^More  to 

• 

«-S 

3.« 

4.« 

3 

«.$ 

t 

t 

«raU  at  tcra 

^tloa 

«f  aapoaai* 

ojet 

O.I 

OJO9 

O.IO 

0.16 

Lipkin  and  Hardy  [11]  measured  the  surface  temperature  of  human  skin 
with  and  without  blood  flow  obstructed  by  exposing  the  surface  to  a 
constant  radiation  of  0.055  cal/sec  cm^ .     Figure  7  shows  the  heating  curves 
of  the  skin  of  the  forearm.     It  is  seen  that  the  temperature  of  the  normal 
skin  increased  to  a  maximum  after  about  60  sec  of  heating  and  then  begin  to 
drop.     However,   the  temperature  of  the  skin  with  blood  flow  obstructed  kept 
increasing.     Thus  the  influence  of  the  blood  flow  which  was  increased  in 
the  normal  skin  due  to  heating  causes  a  decrease  in  skin  temperature  in 
spite  of  continued  irradiation.     This  phenomenon  becomes  apparent  after  the 
first  20  sec.  of  heating. 


Figure  7    Surface  Temperature  of 
Human  Forearm  With  and 
Without  Blood  Flow 
Obstruction  When  Exposed 
to  a  Radiation  of  0.055 


cal/cm^  sec.  From  Lipkin 
and  Hardy  [11] 
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The  effects  of  the  presence  of  edema  fluid  and  the  change  in  blood 
flow  cannot  be  observed  in  Tables  4  and  5.     This  is  because  constant 
thermal  properties  were  assumed  in  the  predictions.     However,  Henriques  and 
Moritz  mentioned  that  even  neglecting  the  variations  of  the  epidermal 
thermal  conductivity  and  other  factors  results  in  uncertains  of  less  than 
0.2*C  in  the  steady  state  temperature  of  the  basal  epidermal  layer. 
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4.       THE  THERMAL  PROPERTIES  OF  SKIN 

The  reaction  of  skin  to  heat  application  depends  significantly  on  its 
thermal  properties  because  these  properties  affect  directly  the  change  of 
skin  temperature  during  heat  transfer  processes.     The  values  of  the 
properties  such  as  conductivity  (k) ,  capacity  (c) ,  density  (p) ,  thermal 
inertia  (kpc)  and  diffusivity  (a  -  k/pc)  of  cutaneous  and  subcutaneous 
tissues  have  been  determined  by  several  investigators.     Some  of  these 
properties  were  found  to  vary  considerably  with  layer,  site,  moisture 
content,  irradiation  and  temperature.     Excised  skin  also  differs  from 
living  skin  because  of  the  increased  blood  flow  due  to  heating. 

Henriques  and  Moritz  [5]  estimated  the  heat  capacity  and  thermal 
conductivity  of  dissected  cutaneous  and  subcutaneous  porcine  tissues.  They 
reported  that  the  average  conductivity  for  epidermis  was  0,0005,  dermis 
0.0009,  fat  0.0004  and  muscle  0.0011  cal/cm  sec'C.     The  average  heat  capa- 
city for  epidermis  was  0.86,  dermis  0.77,  fat  0.55  and  muscle  0.91  cal/g*C. 
The  density  of  epidermis  was  assumed  to  be  0.8  g/cm^.     By  exposing  the  skin 
surface  of  a  pig  to  a  predetermined  temperature  and  maintaining  this 
temperature  at  various  levels  between  45  and  90*C,  they  found  that  the 
thermal  conductivity  of  dermis  increased  two  to  three -fold  after  10  minute 
exposure  due  to  the  presence  of  edema  fluid.     On  the  other  hand,  evidence 
indicated  that  the  thermal  diffusivity  (a  -  k/pc)  was  probably  not  as 
sensitive  as  its  individual  quantities. 

As  reported  in  [5]  and  [11],  the  thermal  conductivity  of  peripheral 
tissues  of  man,  as  determined  by  Lefevre,  and  the  conductivity  of  fat  and 
muscle  for  horse,  pig,  dog  and  cow,  measured  by  Breuer  are  in  good  agree- 
ment with  those  obtained  by  Henriques  and  Moritz.    However,  the  epidermal 


density  from  Leider  and  Bunche  (reported  in  [12])  is  1=1  to  1,2  g/cm^ , 
which  is  about  40  to  50%  higher  than  that  assumed  by  Henriques  and  Moritz . 
Their  results  with  others  are  summarized  in  Table  7, 

Shortly  after  Henriques  and  Moritz,  Hensel  (reported  in  [13])  used  Eq. 
(8)  in  Section  3.1.2  to  calculate  the  thermal  diffusivity  (a)  for  the 
layers  from  the  surface  down  to  the  layers  of  interest  by  measuring  the 
temperature  wave  propagating  into  the  skin  when  contact  heat  was  applied. 
About  one  year  after  Hensel,  Buettner  [13]  used  Eq.   (10)  to  calculate  the 
product  of  kpc  for  surface  heating  of  the  various  parts  of  human  skin. 


Table  7      Thermal  Conductivity  (k) ,  Capacity  (c) 

and  Density  (p)  of  Cutaneous  and  Subcutaneous  Tissues 


lermal  Pronerties 

Enidermis 

Dermis 

Fat 

Muscle 

Remarks 

k 

0.0005 

0.00088 

0.00038 

0.00110 

(1) 

(cal/cm  sec*C) 

0.0005 

/ 

/ 

/ 

(2) 

/ 

/ 

0.00035 

0.00100 

(3) 

/ 

/ 

0.00049 

0.00092 

(4) 

c 

0.86 

0.77 

0.55 

0.91 

(1) 

(cal/g*C) 

0.8 

/ 

/ 

/ 

(1) 

1.1-1.2 

/ 

/ 

/ 

(5) 

(g/cm3) 

/ 

/ 

0.92 

/ 

(6) 

/ 

/ 

/ 

1.27 

(7) 

(1)  From  Henriques  and  Moritz  [5]  for  porcine  tissues 

(2)  Lefevre's  data  for  human  tissues,  reported  in  [11] 

(3)  Breuer's  data  for  porcine  tissues,  reported  in  [5]  and  [11] 

(4)  Hatfield  and  Pugh's  for  human  and  beef  fat,  and  human  muscle, 
reported  in  [11] 

(5)  Leider  and  Bunche 's  reported  in  [12] 

(6)  Handbook  of  Chemistry  and  Physics,  1949 

(7)  From  Lipkin  and  Hardy  [11]  for  human  muscle 
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At  surface  (x  -  0,  T  -  Tg) ,  Eq.   (10)  becomes: 
2 

kpc  ^  ^    ^    ^  (18) 

ir(T     -  T  r 

Thus,  the  product  of  kpc  can  be  evaluated  if  the  temperature  rise  at  the 
surface,  the  heat  flux  and  the  time  are  known.     The  radiation  levels  used 
in  Buettner's  evaluation  ranged  from  0.055  to  0.58  cal/cm^sec.     The  average 
values  of  the  products  figured  out  by  him  were  0.0013,  0.0010  and  0.0023 
cal^/cm^sec°C^  for  human  forearm,  body  and  finger  tip,  respectively.  To 
calculate  the  conductivity  from  these  products,  Buettner  assumed  the  pc  of 
the  skin  to  be  0.8  cal/cm^ 'C,     He  also  calculated  the  conductivity  for  any 
intermediate  layers  based  on  Hensel's  thermal  diffusivity  for  the  layers 
from  the  surface  down  to  any  layers  concerned.     The  conductivity  calculated 
by  Buettner  and  the  diffusivity  determined  by  Hensel  are  shown  in  Table  8. 


Table  8      Thermal  Diffusivity  (k/pc)  and  Conductivity  (k) 
of  Different  Skin  Layers,  from  Buettner  [13] 

 Hensel's  Data  Buettner's  Data  

Layer  ^^P^  Layer  k 

(mm)  (cm^/sec)  (mm)  (cal/cm  sec'C) 


0  -  p. 26  0.0004  0  -  0.26  0.0003 

0  -  0,45  0.0006  0,26  -  0.45  0.0015 

0  -  0.90  0.0010  0.45  -  0.90  0.0024 

0  -  1.30  0.0013  0.90  -  1.30  0.0032 


Lipkin  and  Hardy  [11]  entitled  the  kpc,  the  'thermal  inertia'  for 
surface  heating,  because,  according  to  Equation  (18),  the  greater  the  kpc 
the  lesser  the  surface  temperature  rise  will  be  for  a  given  intensity  of 
heat  flux.     Equation  (18)  was  used  again  by  Lipkin  and  Hardy  to  estimate 
the  kpc.    They  studied  dry  and  moist  excised  tissues,  and  living  skin  with 
blood  flow  intact  and  blood  flow  obstructed.     In  their  studies,  the  levels 


of  radiation  applied  were  around  0.04  to  0,06  cal/cm^sec.     They  found  that 
loss  of  water  in  human  muscle  and  excised  skin  appreciably  decreased  the 
'thermal  inertia',  and  increase  of  blood  flow  in  intact  skin  decreased  the 
rise  of  surface  temperature,  resulting  in  an  increase  in  'thermal  inertia'. 
As  shown  in  Table  9,  the  value  of  the  'inertia'   for  excised  skin  determined 
by  them  ranged  from  about  0.00055  for  dry  skin  to  0.00075  cal^/cm^°C^  sec 
for  moist  skin.     However,  the  'inertia'  for  living  skin  with  blood  flow 
obstructed  was  about  0.0009  cal^/cm^°C^sec .     For  living  skin  with  blood 
flow  intact,   the  'inertia'  varied  from  0.0009  to  0.0040  cal^/cm^'C^sec 
depending  on  the  state  of  dilation  of  blood  vessels .     In  order  to  set  down 
a  simple  value  of  the  inertia  for  surface  heating  for  the  living  skin, 
Lipkin  and  Hardy  suggested  a  value  of  0.0009  to  0.0010  cal^/cm^'C^  sec  for 
short  periods  (20  sec  or  less)  of  heating  by  non-penetrating  radiation  of 
the  dry  skin  of  the  face,  hands,  and  forearms.     While  other  work  by  Hardy, 
Goodell  and  Wolff  (reported  in  [11])  indicated  a  slightly  larger  value  of 
kpc  (0.00125  calVciii^**c2  sec)  for  heating  the  forehead  and  hand  within  3 
sec.     Further,  Hardy  and  Oppel  (reported  in  [11])  obtained  a  value  of 
0.00120  cal^/cm^'C^  sec  during  the  first  30  sec  of  heating  of  the  face. 
Similarly,  the  thermal  conductivity  can  also  be  calculated  from  these 
products.     Lipkin  and  Hardy  suggested  a  value  of  1.20  g/cm^  for  the  density 
and  0.81  cal/g*C  for  the  thermal  capacity. 


Table  9      Thermal  Inertia  (kpc)  of  Human  Tissues, 
from  Lipkin  and  Hardy  [11] 


Tissues 


kpc 

(cal2/cm^ 


sec'C^) 


Muscle 
Fat 

Excised  skin 

Living  skin  (blood  flow  obstruction) 
Living  skin  (blood  flow  intact) 


0.00056  (dry)  to  0.00113  (moist) 
0.00026 

0.00055  (dry)  to  0.00075  (moist) 
0.0009 


0.0009  to  0.004  depending  on  state 
of  vasodilatation 


In  addition  to  the  dependence  on  layer,  site,  moisture  content  and 
blood  flow  rate,  the  thermal  properties  depend  also  on  irradiation  and 
tissue  temperature.     Stoll  and  Greene  [12]  repeated  the  method  to  find  the 
change  of  thermal  inertia  of  the  skin  with  irradiance.     Sites  on  the  volar 
surface  of  the  forearm  were  exposed  to  irradiances  from  0.1  to  0.4 
cal/cm^sec  until  threshold  pain  was  produced.     The  thermal  inertia  at  pain 
threshold  for  different  irradiances  is  given  in  Table  10.     As  shown,  the 
average  kpc  varied  over  a  considerable  range  and  increased  with  each 
increase  in  level  of  irradiance. 


Table  10 


Thermal  Inertia  (kpc)  of  Skin  at 

Pain  Threshold,  from  Stoll  and  Greene  [12] 


Irradiance 
(cal/cm^sec) 


0.10 
0.15 
0.20 
0.30 
0.40 


96 
110 
127 
139 
159 


75  -  117 

89  -  142 

95  -  158 

115  -  168 

144  -  181 


Finally,  Weaver  and  Stoll  [10]  studied  the  variations  of  the  thermal 
conductivity  with  tissue  temperature  for  both  heating  and  cooling  phases 
using  Eq.    (11).     They  assumed  that  pc  =  1  cal/cm^°C  for  human  skin  so  that 
Q  =:  k  and  the  conductivity  remains  constant  from  the  surface  inward  to  the 
basal  layer,  the  interface  between  the  dermis  and  the  epidermis.     Thus,  Eq . 
(11)  can  be  solved  for  either  the  thermal  diffusivity  or  the  conductivity 
for  any  given  experimental  time -temperature  history  at  80  fxm,   the  best 
value  for  the  depth  of  the  basal  layer  in  the  volar  surface  of  the  forearm 
([5,  10]).     Some  of  their  results  are  given  in  Fig.  8  in  which  each  of  the 
curves  can  be  separated  into  two  curves.     The  upper  curve  sloping  downward 
to  the  right  represents  a  continuous  heating  to  a  peak  temperature  of  about 
52  to  57"C  depending  on  the  radiation  levels.     The  lower  curve  sloping 
downward  to  the  left  represents  a  continuous  cooling  from  the  peak  tempera - 
ure  back  to  44*C.     The  exposure  time  for  each  curve  is  the  shortest  time 
required  to  produce  blistering  within  24  hours  (blistering  may  be  delayed 
in  some  instances,  see  Section  5.1),  and  ranges  from  5.56  to  33.80  sec 
depending  on  the  radiation.     As  seen  in  the  figure,  the  thermal  conduc- 
tivity decreases  continuously  when  the  tissue  temperature  increases  during 
the  heating  phase.     During  the  cooling  phase,  the  conductivity  decreases 
with  the  tissue  temperature.     The  greatest  variations  occur  at  temperatures 
below  44*C,  but  little  variation  occurs  at  temperatures  above  44"C,  The 
effects  of  the  presence  of  edema  fluid  and  the  increase  of  blood  flow  due 
to  continuous  heating  were  not  mentioned  in  their  paper. 

Another  thermal  property  that  has  not  yet  been  discussed  is  the 
absorption  coefficient  (u) .     This  coefficient  is  significant  only  when 
penetrating  radiation  is  considered.     Buettner  ([6],   [13])  found  that  the 


Figure  8    Variations  of  the  Thermal  Conductivity  of  Human  Skin  with  Tissue 
Temperature  at  80  /xm  deep.     From  Weaver  and  Stoll  [10] 

absorption  coefficient  of  human  skin  was  between  100  and  1000  cm'^  for  the 
spectrum  of  infrared  radiation.     Since  this  coefficient  is  very  large, 
infrared  radiation  does  not  penetrate  the  skin.     The  solar  rays,  however, 
showed  a  i/  near  10  cm"^.    According  to  Refs.  6  and  14,  the  spectrvim  of  the 
radiation  from  an  atomic  bomb  blast  may  be  considered  close  to  that  of  the 
sun.     Thus  a  value  »/  -  10  cm'^  can  be  assumed  for  this  case. 
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5.       THE  QUANTITATIVE,  QUALITATIVE  AND  SENSATIONAL  EFFECTS  OF  HEAT  ON  THE 
SKIN 

Reactions  of  skin  to  heat  applications,  ultimately  result  in  cutaneous 
burns.     The  severity  of  burns  may  be  characterized  quantitatively  according 
to  the  depth  to  which  the  tissue  has  been  destroyed,  or  qualitatively 
according  to  the  nature  of  the  changes  that  have  occurred,  or  even 
sensationally  according  to  the  intensity  of  pain  induced. 

5 . 1    General  Consideration  of  the  Quantitative  and  Qualitative  Effects 

Moritz  [15]  studied  both  the  quantitative  and  qualitative  effects  of 
heat  on  the  skin.     He  indicated  that  the  mildest  form  of  epidermal  injury 
produced  by  burning  was  latent  in  the  sense  that  alteration  in  cell 
structure  were  unable  to  be  recognized,  although  temporary  impairment  in 
the  adhesion  between  epidermis  and  dermis  were  observed.     Such  injuries 
were  reversible  and  the  time  required  for  their  recovery  increased  with  the 
time  required  for  their  production. 

The  first  evidence  of  irreversible  thermal  injury  of  the  epidermis  was 
a  redistribution  of  water  and  solids  within  the  nuclei  of  the  epidermal 
cells.    As  a  result,  the  nuclei  swelled,  first  in  the  intermediate  and 
later  in  the  deepest  layer.     Further  injury  resulted  in  swelling  and 
disintegration  of  the  cytoplasm  (a  complex,  jelly-like  colloidal  substance 
constituting  the  living  matter  outside  a  cell  nucleus)  of  the  basal  cells 
and  pyknosis  (thickening  and  shrinking  through  degeneration)  of  nuclei 
throughout  the  entire  thickness  of  the  epidermis.    When  the  temperature 
within  the  epidermis  was  less  than  about  55*C,  further  primary  thermal 
alterations  of  the  epidermal  cells  were  negligible.    When  the  temperature 
exceeds  approximately  55 'C,  coagulation  occurred.     Coagulative  changes 


29 

might  mask  both  nuclear  swelling  and  cytoplasmic  disintegration.  Further 
increases  in  temperature  resulted  in  desiccation  and  eventually  in 
carbonization. 

The  earliest  visible  alternative  change  in  the  dermis  was  constriction 
of  the  superficial  blood  vessels.     Normally,  this  was  followed  immediately 
by  dilation.     Thermal  dilation  of  superficial  blood  vessels  lead  to 
increased  vascular  permeability  and  edema  formation.     When  a  heat  exposure 
was  sufficient  to  damage  the  attachment  between  epidermis  and  dermis ,  the 
amount  of  edema  fluid  that  collected  between  them  was  usually  sufficient  to 
elevate  and  to  form  a  blister.     Therefore,  blistering  customarily  indicated 
destruction  of  epidermis.     However,  absence  of  blisters  did  not  always 
indicate  epidermal  survival  because  there  were  several  instances  in  which 
blistering  was  either  delayed  or  prevented. 

Delayed  blistering  was  most  often  observed  after  long  exposures  at  low 
temperatures  or  after  flash  exposures  at  high  temperatures .     In  both  cases , 
the  exposures  were  either  too  low  in  intensity  or  too  short  in  duration  to 
cause  sufficient  vascular  damage.    As  a  result,  long  hours  (usually  within 
24  hours)  were  required  to  collect  enough  fluid  under  the  damaged  epidermis 
to  form  a  blister.     Also,  after  extreme  exposures,  the  dermis  and  its 
superficial  blood  vessels  were  almost  immediately  coagulated.     In  this 
case,  blood  flow  through  superficial  capillaries  was  stopped  and  edema 
developed  was  too  deep  to  produce  blistering. 

Like  the  epidermis,  coagulative  changes  occurred  in  the  dermis  at  a 
temperature  of  approximately  55 or  above.     Further  increases  in  dermal 
temperature  might  result  in  desiccation  and  carbonization.     Thus,  for  a 
exposure  of  insufficient  intensity  to  raise  the  epidermal  or  dermal 


temperature  above  55**C,  coagulation  did  not  occur  even  though  deep 
destruction  might  be  produced  if  sufficiently  prolonged,     while  a 
relatively  shallow  bum  might  show  coagulation  and  even  carbonization  if 
the  exposure  had  been  intense  and  short. 

The  depth  to  which  the  tissue  was  destroyed  (the  quantitative 
reaction) ,  on  the  other  hand,   increased  with  both  the  duration  and 
intensity  of  exposures.     But  the  duration  and  intensity  of  exposures  to 
produce  similar  quantitative  reactions  were  inversely  related.     In  this 
case,   the  quantitative  results  of  a  short  exposure  of  high  intensity  may  be 
similar  to  those  of  a  long  exposure  at  low  intensity.     As  discussed 
previously,  a  high  intensity  exposure  usually  results  in  a  coagulative  type 
of  bum  and  a  low  intensity  results  in  a  non- coagulative  type  of  burn. 
Apparently,  quantitative  similarity  between  reactions  does  not  necessarily 
mean  that  the  reactions  are  qualitatively  similar.     Therefore,  the 
superficial  appearance  of  a  bum  is  not  a  reliable  criterion  for  judging 
the  depth  to  which  the  tissue  may  have  been  destroyed  [15]. 

In  order  to  classify  the  severity  of  thermal  injury  quantitatively, 
cutaneous  bums,  according  to  [15,  16]  were  characterized  as  first,  second, 
or  third  degree.     First  degree  bums  referred  to  the  cutaneous  reactions 
that  indicated  incomplete  destruction  of  the  epidermis.     Second  degree 
bums  were  defined  as  the  reactions  that  indicated  complete  destruction  of 
the  epidermis  without  a  significant  degree  of  injury  to  the  dermis. 
Reactions  indicative  of  a  complete  destruction  of  the  epidermis  and  a 
serious  damage  to  the  dermis  were  designated  as  third  degree. 
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5 . 2    General  Consideration  of  the  Sensational  Effect 

Pain  is  the  primary  sensational  effect  of  heat  on  skin.     The  intensity 
of  pain  caused  by  cutaneous  burning  can  be  explained  according  to  the 
concept  introduced  by  Weber  (reported  in  [17]).     The  Weber  Law  states  that 
as  the  intensity  of  a  stimulus  increases,  the  ability  to  discriminate 
differences  in  intensity  of  stimulus  decreases  proportionately.     In  other 
words,  as  the  intensity  (I)  of  a  stimulus  increases,   the  minimum  change  in 
intensity  (Aim) »  which  can  just  barely  be  distinguished  or  felt  when  added 
to  or  subtracted  from  I,  increases  proportionately.     That  is: 


in  which  AI^/I  is  known  as  the  Weber  ratio  and  Ci  is  a  constant. 

Fechner  extended  Equation  (19)  to  relate  the  sensation  to  the 
stimulation  by  assuming  that  for  a  given  intensity  of  stimulus  (I) ,  the 
factor  by  which  the  change  in  intensity  of  stimulus  (AI)  exceeds  the  "just 
noticeable  difference"  in  intensity  (AI^)  is  proportional  to  the  change  in 
sensation  (AS).     That  is: 


where  k  -  C1C2.  Upon  integration  from  the  threshold  stimulus  intensity  Iq 
to  an  intensity  I,  the  corresponding  intensity  of  sensation  S  is: 


(20) 


Substituting  Equation  (20)  into  Equation  (19)  gives: 


J-  -  kAS 


(21) 


S  -  So  +  k  log  I/Io 


(22) 


where  Sq  -  threshold  sensation. 
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Thus,  Equation  (22)  can  be  used  to  evaluate  the  intensity  of  pain  sensation 
for  any  types  of  thermal  stimulus , 

A  scale  of  pain  intensity  has  been  established  by  Hardy,  Wolff  and 
Goodell  [17].     They  measured  the  "just  noticeable  difference"  in  intensity 
of  thermal  radiation  stimulus  (Aim)  each  increase  in  radiation  level 

starting  at  the  intensity  at  which  threshold  pain  occurred  to  the  intensity 
at  which  further  increase  in  radiation  caused  no  perceptibly  greater  pain. 
The  pain  induced  by  this  maximal  magnitude  of  radiation  is  therefore  a 
"ceiling"  pain.     They  found  that  there  were  21  discriminate  intensities  of 
pain  between  the  pain  threshold  and  the  ceiling  pain.     Hardy  et  al.,  then, 
set  up  a  pain  intensity  scale  in  terms  of  just  noticeable  difference  (JND) 
according  to  the  assumption  of  Fechner  that  the  intensity  of  sensation 
evoked  by  a  given  stimulus  was  equivalent  to  the  number  of  the 
discriminable  steps  from  the  threshold  sensation.     Thus,  a  value  of  0  JND 
indicates  the  beginning  of  the  scale  or  no  pain  and  a  value  of  21  JND 
indicates  the  maximum  of  pain  intensity.     Another  unit  of  pain  intensity 
called  "dol"  was  also  proposed  by  them.     They  defined  1  dol  <-  2  JND. 
Figure  9  indicates  the  pain  intensity,  in  terms  of  both  dol  and  JND, 
induced  by  heating  the  forehead  and  forearm  for  a  duration  of  3  sec.  As 
shown,  a  pain  threshold  was  produced  at  about  0.22  cal/cm^sec  for  a  3  sec 
exposure.    When  the  intensity  was  increased  above  0.48  cal/cm^sec,  it  was 
not  until  the  stimulus  had  been  increased  to  0.68  cal/cm^sec  that  a  barely 
perceptible  difference  from  0.48  cal/cm^sec  could  be  noted.  Although 
intensity  beyond  0.68  cal/cm^sec  could  be  distinguished  from  0.48 
cal/cm^sec,  it  was  not  possible  to  distinguish  any  stimuli  greater  than 
0.68  from  0.68  cal/cm^sec.    Therefore,  the  sensation  evoked  by  a  stimulus 
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Figure  9    Pain  Intensity  Induced  by  Heating  the  Forehead  and  Forearm  for  a 
Duration  of  3  Seconds.     From  Hardy,  Wolff  and  Goodell  [17] 

of  about  0.68  cal/cm^sec  was  the  ceiling  pain  and  no  greater  pain  could  be 
felt  beyond  this  point.     Hardy  et  al.  also  indicated  that  beginning  at  the 
threshold  the  type  of  cutaneous  pain  was  a  pricking  sensation  and  after  the 
stimulus  intensity  was  increased  to  about  0.32  cal/cm^sec,  a  burning 
quality  of  pain  was  added  to  the  total  sensation. 

5 . 3    Thresholds  of  Bums.  Blistering  and  Pain  as  Related  to  the  Intensity 
and  Duration  of  Thermal  Exposure 

The  time -temperature  thresholds  for  first  and  second  degree  bums  in 
pigs  and  men  have  been  studied  by  Moritz  and  Henriques  [16]   (pigs  were  used 
because  porcine  skin  is  similar  to  human  skin) .    Although  it  was  found  that 
a  pig  did  not  sweat  like  a  man,  and  that  there  were  certain  qualitative 
differences  in  the  responses  of  human  and  porcine  skin  to  heat  application, 


there  were  no  significant  quantitative  differences  in  their  sensitiveness 
to  thermal  injury.     Moritz  and  Henriques  also  discovered  that  porcine  skin 
was  not  uniformly  sensitive  to  thermal  injury.     The  skin  of  the  neck  and 
the  mid-portion  of  the  back  was  found  to  be  less  sensitive  than  that  of  the 
lateral  portion  of  the  shoulders,  thorax  and  abdomen.     While  the  most 
sensitive  areas  were  the  portions  on  the  ears,   thighs,  buttocks  and  ventral 
surface.     For  human  skins,  Buettner  [13]  illustrated  that  the  forearm 
appeared  to  be  less  susceptible  to  thermal  injury  than  any  other  parts  of 
the  body  (e.g.  forehead,  chest,  upper  arm,  palm,  leg,  etc.). 

Moritz  and  Henriques  obtained  the  overall  results  of  the  time- 
temperature  thresholds  at  which  cutaneous  burns  occurred,  resulting  from 
179  exposures  of  porcine  skin  and  33  exposures  of  human  skin  on  different 
parts  of  the  body.     Their  results  are  plotted  in  Figure  10  in  which  the 
surface  temperature  was  attained  by  contacting  the  skin  surface  directly 
with  a  rapidly  flowing  stream  of  hot  water  or  oil  maintained  at  constant 
temperatures.     The  crosses  fitted  by  the  solid  line  represent  the  minimum 
exposures  that  caused  complete  destruction  of  porcine  epidermis  (second 
degree  threshold) .    Any  exposures  situated  above  the  solid  line  are  supra- 
threshold  and  those  sitxiated  below  are  sub- threshold.     The  broken  line 
represents  the  maximum  exposures  that  failed  to  destroy  porcine  epidermis 
(first  degree  threshold).    Thus  exposures  lying  below  the  broken  line  cause 
no  appreciable  injury.    Reactions  of  human  skin  are  indicated  by  the 
circles.     The  solid  circles  indicate  the  second  degree  threshold,  while 
the  open  circles  represent  the  first  degree  threshold.    As  seen,  there  are 
no  significant  differences  in  the  time- temperature  thresholds  between  human 
and  porcine  skins. 


Figure  10  indicates  an  inverse  relationship  between  the  surface 
temperature  and  the  time  required  to  produce  first  or  second  degree  burns. 
The  minimum  surface  temperature  required  to  produce  a  threshold  burn  is 
about  UWC  and  the  time  is  about  6  hours.     At  surface  temperatures  between 
UU^C  and  51°C,  the  rate  at  which  burning  occurs  increases  logarithmically 
with  a  linear  increase  in  skin  temperature.     Moritz  and  Henriques  mentioned 
that  in  the  case  of  surface  temperature  under  Sl^C  the  total  exposure  time 
required  to  destroy  the  epidermis  was  (a)  so  long  that  the  time  required  to 
bring  the  temperature  at  the  basal  layer  to  a  steady  state  was 
comparatively  negligible,  and  (b)  essentially  equal  to  the  total  duration 
of  the  steady  thermal  state  within  the  epidermis. 
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Figure  10      Time -Surface  Temperature  Thresholds  at  Which  Cutaneous 
Burning  Occurs.     From  Moritz  and  Henriques  [16]. 


At  surface  temperatures  above  Sl'C,   the  rate  of  injury  increases 
rapidly  and  the  pre-burn  time  is  too  short  to  allow  the  epidermis  to  reach 
steady  state.     Consequently,   the  relationship  between  the  surface 
temperature  and  the  pre-burn  time  is  no  longer  linear.     Above  70°C,  a 
second  degree  burn  can  be  produced  in  less  than  1  second. 

The  same  type  of  temperature  -  time  threshold  curve  as  presented  by 
Moritz  and  Henriques  has  been  confirmed  by  other  investigators.     They  found 
that  all  these  data  could  be  represented  by  a  simple  exponential  function 
[18]: 

(T^ef  ■  ^> 

t  -  t^^,e  (23) 

where  t        and  T     _  are  the  reference  time  and  temperature  respectively, 
ref  ref  ^  v  j 

Similarly,   the  thresholds  of  pain  and  blistering  produced  by  exposures 
of  thermal  radiation  have  also  be  determined  [12,   19,   20].   Figure  11  illus- 
trates the  tolerance  times  for  different  levels  of  thermal  energy  which  are 
plotted  on  log  log  coordinates  [20].     In  this  figure,   the  tolerance  time 
represents  the  time  required  to  produce  threshold  pain  and  minimal 
blistering  (i.e.   several  small  blisters  appearing  within  24  hours).  These 
results  were  obtained  from  exposures  of  the  volar  surface  of  human  forearm 
with  an  initial  skin  temperature  of  32.5°C.     As  seen,  a  similar  inverse 
relationship  between  the  tolerance  time  and  the  intensity  of  absorbed 
energy  appears  in  both  the  threshold  curves  of  pain  and    blistering.  The 
area  between  the  2  curves  represents  the  region  of  suprathreshold  pain  and 
reversible  tissue  damage.     The  area  above  the  blister  curve  represents  the 
region  of  irreversible  injury  and  below  the  pain  curve,  the  safety  zone. 
The  survival  curve  in  between  the  pain  and  blister  curves  is  of  no  interest 
at  the  present  time  and,  so,  will  not  be  discussed. 
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Figure  11      Human  Skin  Tolerance  Time  to  Absorbed  Thermal  Energy. 
From  Stoll  and  Chianta  [20] 

As  shown  in  Figure  11,  the  intensity  and  duration  at  which  threshold 

blistering  occurs  range  from  0.1  cal/cm^  sec  in  about  33  sec  to  1.3 

cal/cm^  sec  in  less  than  1  sec.     For  the  pain  threshold,  the  range  is 

between  0.1  cal/cm^  sec  in  approximately  13  sec  and  0.64  cal/cm^  sec  in 

about  1  sec. 
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Although  threshold  pain  and  blistering  are  functions  of  both  exposure 
time  and  radiation  level,  threshold  pain  occurs  at  some  definite  tempera- 
ture, almost  independently  of  the  time  required  to  attain  this  temperature, 
and  also  independently  of  the  site  of  the  body,  whereas  blistering  depends 
strongly  on  both  temperature  and  time  [12,   13,  19]. 

As  reported  by  Hardy,  Jacobs  and  Meixner  [19],   threshold  pain  was 
elicited  at  a  surface  temperature  of  about  45.7°C  (Table  11)  and  ceiling 
pain  at  about  60*C  (Figure  12)  regardless  of  the  time  of  exposure  and 
intensity  of  stimulus.     In  this  case,   the  intensity  of  pain  can  be 
expressed  simply  in  terms  of  the  surface  temperature  as  shown  in  Figure  12 . 

However,  Stoll  and  Greene  [12]  illustrated  that  the  skin  surface 
temperature  at  pain  threshold,  instead  of  being  a  constant  as  shown  in 
Table  11,  tended  to  slightly  increase  at  higher  radiation  levels  due  to  the 
greater  temperature  gradient  within  the  dermis  at  the  higher  intensities, 
but  the  threshold  temperature  of  the  pain  receptors  at  an  average  depth  of 
about  0.2  mm  was  invariable  (about  43. 2 "C)  throughout  the  range  of 
irradiances  between  0.1  and  0.4  cal/cm^  sec.     Since  the  0.2  mm  depth  was 
calculated  on  the  basis  of  the  data  obtained  from  exposures  of  the  forearm, 
Stoll  and  Greene  pointed  out  that  (a)  this  value  was  much  too  thick  for 
epidermis  on  the  forearm,  the  effective  pain  receptors  were  therefore 
located  in  the  dermis,  and  (b)  the  significance  of  this  location  was  that 
it  was  theoretically  feasible  to  destroy  completely  the  epidermis  without 
causing  any  pain,  as  long  as  the  temperature  gradient  through  the  skin  was 
such  that  the  epidermal  temperature  was  maintained  at  a  threshold  or 
suprathreshold  level  while  the  pain  receptors  were  maintained  at  a  sub- 
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Table  11      Stimulus  Intensity,  Exposure  Time  and  Skin  Surface  Temperature 
at  Pain  Threshold.     From  Hardy,  Jacobs  and  Meixner  [19] 
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Figure  12      Relationship  Between  Pain  Intensity  and  Skin  Surface 

Temperature  Obtained  by  Heating  the  Surface  of  the  Forehead. 
From  Hardy,  Jacobs  and  Meixner  [19]. 


threshold  level.  This  is  probably  the  reason  why,  as  mentioned  in  [19], 
thermal  injury  without  pain  and  pain  without  apparent  injury  are  common 
experiences , 

Before  Stoll  and  Greene,  Buettner  [13]  indicated  that  'unbearable' 
pain,  i.e.  not  the  least  perceptible  pain  as  discussed  previously,  occurred 
at  a  threshold  temperature  of  about  44,8''C  at  a  depth  of  0.1  mm  rather  than 
0.2  mm.     While  other  investigators  argued  that  free  nerve  endings  in  the 
epidermis  might  mediate  pain  [12]. 

5.4    Factors  Affectine  the  Effects  of  Thermal  Exposure  on  the  Skin 

The  possiblity  of  altering  the  vulnerability  of  skin  to  thermal  injury 
by  changing  the  initial  skin  temperature,  blood  pressure  and  water  content 
of  the  skin  has  also  been  studied.     It  was  found  that  a  change  in  initial 
skin  temperature  did  not  alter  the  pain  threshold  temperature  but  the  times 
required  to  produce  threshold  pain  and  burns  were  significantly  affected 
during  radiation  exposures  [12,  13,  19].     Stoll  and  Greene  [12]  demon- 
strated that  the  rate  at  which  burning  occurred  were  markedly  increased 
even  though  the  initial  temperature  was  increased  by  only  l^C  for  the  same 
intensity  of  irradiance.     The  initial  skin  temperature  measured  by  them  was 
32.5'C.     Buettner  [13]  indicated  that  precooling  the  skin  from  an  average 
initial  temperature  of  about  33*C  significantly  delayed  the  production  of 
pain  and  bums.     In  the  case  of  direct  thermal  contact,  Henriques  [21] 
mentioned  that  the  effect  of  initial  temperature  on  the  rate  of  change  of 
basal  layer  temperature  was  relatively  insignificant.     He  suggested  a  value 
of  35°C  for  the  initial  skin  temperature. 
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Changes  in  blood  pressure  by  compression  of  the  skin,  on  the  other 
hand,  showed  little  effect  on  the  susceptibility  to  thermal  injury  although 
it  altered  the  surface  color  of  the  burns  [15,  16].     Alteration  of  water 
content  by  wetting  the  skin,  according  to  [13],  exerted  no  appreciably 
influence  on  the  threshold  time- intensity  relationships  at  intensities  of 
radiation  above  0.25  cal/cm^  sec.     At  about  0.055  cal/cm^  sec,  wet  skin  was 
found  to  be  less  sensitive  than  dry  skin  nevertheless  the  threshold  pain 
temperature  remained  unchanged. 

Another  factor  which  is  worthy  of  mention  is  the  cumulative  effects  of 
repeated  subthreshold  exposures.     Although  a  single  subthreshold  exposure 
causes  no  recognizable  epidermal  injury,  latent  injury  may  be  sustained. 
When  several  subthreshold  exposures  are  repeatedly  applied  to  the  same 
area,  latent  injury  produced  by  each  exposure  can  be  accumulated  to  become 
recognizable  injury  if  the  intervals  between  exposures  are  shorter  than  the 
time  required  for  their  recovery. 

Table  12  demonstrates  this  effect  on  the  porcine  skin.     Exposures  were 
made  with  a  running  stream  of  hot  water  at  49 "C  and  at  atmospheric  pressure 
[16].     Exposures  identified  by  Reference  Nos.  1  to  18  are  single  exposures 
for  control  purposes.    As  seen,  single  exposures  for  less  than  7  minutes 
cause  no  recognizable  epidermal  injury.    Repeated  exposures  are  Indicated 
by  Nos.  19  to  39.    Though  a  single  3  minute  exposure  causes  no  recognizable 
injury,  three  such  exposures  each  separated  by  recovery  periods  less  than 
48  minutes  have  the  same  effect  as  a  single  continuous  9  minute  exposure. 
For  intervals  between  exposures  greater  than  48  minutes,  appreciable 
recovery  was  observed.     Complete  recovery  requires  2  to  4  hours  as 
indicated  by  Experiments  30  and  31. 
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Table  12      The  Cumulative  Effects  of  Repeated  Subthreshold  Thermal 

Exposures*  on  Porcine  Skin.     From  Moritz  and  Henriques  [16] 
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*An  exposures  were  made  to  water  at  49°C 

**The  sum  of  the  morphological  changes  indicative  of  cell 
death  and  caused  by  the  progressive  degradative  action 
of  enzymes 


Comparing  the  exposures  from  Nos.  34  to  39,  one  finds  that  the  time 
required  for  recovery  from  latent  injury,  as  would  be  expected,  increased 
with  the  time  of  exposure.     Therefore,  the  degree  of  recovery  depends 
partly  on  the  length  of  the  recovery  period  and  partly  on  the  duration  of 
exposure  which  causes  unrecognizable  injury. 


6.       PREDICTION  OF  CUMULATIVE  TISSUE  DAMAGE  INCURRED  DURING  A  BURN  AND  ITS 
RELATION  TO  PAIN  INTENSITY 


Before  predicting  the  tissue  damage  and  its  relation  to  pain  inten- 
sity,  it  is  necessary  to  understand  briefly  the  physical  and  chemical 
processes  that  take  part  in  cell  death  during  a  burn.     These  processes, 
according  to  Henriques  [21] ,  can  be  basically  classified  into  three  general 
groups : 

(a)  Thermal  denaturation  of  proteins  -  Living  cells  are  primarily  protein 
in  nature.     These  proteins  which  participate  in  the  maintenance  of 
normal  cell  life  are  very  sensitive  to  heat. 

(b)  Other  changes  in  metabolic  processes  regardless  of  their  effect  on 
protein  -  For  instance,  increase  in  tissue  temperature  may  alter  the 
rate  of  diffusion,   formation  and  degradation  of  chemical  reactants. 
This  will  lead  to  a  change  in  the  concentrations  of  some  of  the 
individual  constituents.     Consequently,  the  entire  metabolic 
equilibrium  may  be  disturbed. 

(c)  Alterations  in  the  physical  characteristics  of  cells  -  This  group 
includes  all  physical  phenomena  that  are  characteristic  of  protoplasm 
(a  complex,  jellylike  celloidal  substance  constituting  the  living 
matter  of  a  cell)  but  are  not  essentially  affected  by  the  thermal 
degradation  of  proteins,  for  example,  diffusion  through  liquids  and 
membranes,  changes  in  viscosity,  rigidity  and  liquefaction.     All  these 
biophysical  processes  are  potentially  able  to  cause  cell  death. 

For  any  of  the  above  processes  to  take  place,  the  constituents 
involved  should  possess  an  energy  content  which  is  sufficiently  high  to 
activate  the  process.     The  minimum  amount  of  energy  required  for  the 
process  to  proceed  is  known  as  the  activation  energy.     The  rate  of  the 
process  depends  directly  on  the  fraction  of  the  constituents  which  have  an 
energy  content  at  least  equal  to  the  activation  energy.     In  most  chemical 
and  physical  processes,  this  fraction  (f)  based  on  the  Maxwell -Boltzman 


energy  distribution  law  [22]  is: 
^  ^  ^-AE/RT 


(24) 


where  AE  -  activation  energy,  cal/mole 
R    -  gas  constant  *»  2  cal/mol  K 
T    -  temperature  in  K 

Henriques  assumed  that  the  rate  of  any  conceivable  process  that  might 

result  in  cell  death  could  be  expressed  in  a  formula  similar  to  that  for 

most  chemical  rate  processes,  the  Arrhenius  equation  [22).  Thus: 

dfl  _  p  ^-  AE/RT  (25) 
dt 

where  dQ/dt  -  rate  of  cell  damage,  sec"^ 

Q  -  an  arbitrary  function  to  represent  the  amount  of  cell  damage 
T  -  tissue  temperature  in  K 
P  -  constant,  sec 

Then  the  cumulative  damage  is: 

Q  -  P  ;t  e-^E/RT  ^26) 
As  described  in  Section  5.3,  when  the  skin  surface  is  maintained  at  a 
temperature  less  than  50*C  during  exposures  by  direct  thermal  contact,  the 
total  exposure  time  required  to  reach  second  degree  burn  threshold  is 
essentially  equal  to  the  total  duration  of  the  steady  thermal  state  within 
the  epidermis.     In  this  case,  the  temperature  at  the  basal  layer  can  be 
considered  invariant  and  equal  to  the  steady  state  temperature.  Equation 
(26)  can  be  integrated  to  obtain: 

0  -  P  t  e-^^""  (27) 
or 

In  t  -  In  0  -  In  P  +  AE/RT  (28) 
Henriques  assigned  O  -  1  to  represent  the  second  degree  burn  threshold 
(damage  just  sufficient  to  cause  complete  destruction  of  the  epidermis). 
Together  with  the  data  given  in  Figure  10  for  surface  temperatures  less 
than  50*C,  the  constant  P  and  the  overall  activation  energy  AE  including 
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all  the  processes  involved  in  cell  damage  can  be  obtained  from  the 
intercept  and  slope  of  the  plot  In  t  vs.  1/RT  (Equation  (28)).  These 
values  were  found  to  be : 

P  -  3.1  X  lO^S  sec'^  (29) 
AE  -  150,000  cal/mole  (30) 

Substituting  Equations  (29)  and  (30),  and  the  data  pertaining  to  the 
first  degree  burn  threshold  in  Figure  10  back  to  Equation  (28) ,  the  best 
fitted  damage  function  representing  the  maximum  damage  condition  before 
occurrence  of  irreversible  epidermal  injury  is: 

Q  «  0.53  (31) 

Although  the  values  of  P,  AE,  and  0  were  obtained  by  using  Equation 
(27)  and  a  limited  portion  of  the  experimental  data,  Henriques  found  that 
Equation  (26),  with  these  values  of  P  and  AE  provided,  could  accurately 
predict  the  degree  of  epidermal  injury  in  terms  of  the  above  defined  Q  over 
the  entire  experimental  temperature  range  for  both  direct  contact  and 
ambient  hot  air  exposures,  as  long  as  the  temperature  history  at  the  basal 
layer  was  known. 

Henriques  also  found  that  of  all  of  the  physical  and  chemical 
processes  potentially  capable  of  causing  cell  damage,  only  the  process  of 
protein  denaturation  had  an  activation  energy  roughly  equal  to  150,000 
cal/mole.    The  other  processes  had  far  lower  activation  energy.     In  this 
case,  Henriques  concluded  that  thermally  induced  epidermal  injury  was 
mainly  the  result  of  thermal  alterations  of  those  cellular  proteins  whose 
rate  of  denaturation  corresponded  to  an  activation  energy  of  at  least 
150,000  cal/mole. 


Stoll  and  Greene  [12]  used  the  rate  process  model  as  derived  by 
Henriques  to  estimate  skin  damage  due  to  thermal  radiation.     They  measured 
the  temperature  history  at  the  skin  surface  for  each  irradiance  used. 
Based  on  these  temperatures,  the  tissue  temperature  at  a  depth  of  80  /xm, 
approximately  the  basal  layer,  at  every  instant  in  time  was  calculated 
using  Equation  (8)  in  Section  3.1,2.     Although  Equation  (8)   is  not  for 
thermal  radiation  exposure,  Stoll  and  Greene  claimed  that  this  equation  was 
equally  valid  when  the  surface  temperature  at  every  instant  was  known.  The 
calculated  tissue  temperature  was  then  substituted  into  Equation  (25)  to 
solve  for  the  rate  of  damage  dfl/dt.     By  plotting  dfl/dt  against  the  time  at 
which  the  tissue  temperature  occurred,  the  area  under  the  resultant  curve 
was  then  the  cuunulative  damage  of  the  epidermis.     For  the  production  of 
threshold  blistering  which  is  equivalent  to  the  second  degree  burn 
threshold,  Stoll  and  Greene  found  that  the  total  damage,  including  the 
damage  done  during  both  heating  and  cooling  phases,  was  about  two- tenths  of 
the  value  obtained  by  Henriques  for  the  same  injury. 

Upon  analysis  of  all  the  possible  reasons  for  this  discrepancy,  Stoll 
and  Greene  commented  that  since  the  steady  state  was  never  achieved  at  skin 
temperatures  above  50*C  in  conductive  heating,  there  might  be  a  discontin- 
uity in  the  damage  rate  function  which  was  inconspicuous  by  the  transition 
from  the  steady  state  to  the  \insteady  state  conditions.     They  further 
cooomented  that  the  damage  imposed  after  the  exposure  was  removed  should  not 
be  overlooked,  especially  at  high  intensity.  As  they  mentioned,  about  35% 
of  the  total  damage  was  incurred  during  the  return  of  the  tissue  tempera- 
ture to  the  initial  level  at  intensity  of  0.4  cal/cm^  sec  (about  10%  at  0.1 
cal/cm^  sec) .     In  order  to  fit  all  their  data  to  the  rate  process  model  to 


yield      -  1  to  represent  the  total  damage  performed  during  heating  and 
cooling  at  the  blister  threshold,  they  obtained  the  curve  shown  in  Figure 
13.  which  indicates  the  log  of  the  damage  rate  against  the  tissue 
temperature  at  80  /xm.     The  damage  rate  at  tissue  temperature  for  conduct! 
heating  from  Henriques  is  also  presented. 


44       ««        4«        so        ftZ        »4        »«  M 

TEMPCRATUMC  °C 

Figure  13      The  Damage  Rate  Against  the  Tissue  Temperature  at  80  /xm  Such 

That  0  -  1  at  Threshold  Blistering.  From  Stoll  and  Greene  [12]. 

To  relate  pain  to  tissue  damage,  Hardy  et  al.   [19]  analyzed  the  basic 
event  underlying  'noxious  stimulation'  resulting  from  thermal  impact 
('noxious  stimulation'  refers  to  the  effect  evoked  at  the  pain  fiber 
endings,  i.e.  pain,  reflex  response,  etc.,  rather  than  to  the  type  of 
stimulation  used  to  produce  these  responses).     They  discovered  that  tissue 
damage,  protein  denaturation,  pain  and  reflex  responses  (e.g.  withdrawal  of 
exposed  area  in  man;  contraction  of  muscle  in  guinea  pig;  and  tail  flick  in 


rat  in  response  to  thermal  stimulus)  all  began  at  a  threshold  temperature 
of  about  AS'C.     This  temperature  held  in  spite  of  the  type,   intensity  and 
duration  of  thermal  stimulus  and  also  regardless  of  the  initial  skin 
temperature.     Hardy  et  al.   inferred  that  the  same  process  underlined 
noxious  stimulation  and  tissue  damage,   that  is  thermal  denaturation  of 
cellular  proteins.     Although  thermally  induced  pain  and  damage  can  be 
explained  as  the  results  of  the  same  process,  pain  depends  only  on  the 
temperature  to  which  the  skin  is  raised  while  tissue  damage  depends  on  both 
temperature  and  time.     From  this  standpoint,  Hardy  et  al .   further  inferred 
that  noxious  stimulation  of  the  skin  eliciting  pain  was  dependent  on  the 
rate  of  damage  (dft/dt)  rather  than  the  amount  of  damage  (Q)   induced  during 
a  burn.     This  is  because  Equation  (25)  is  a  function  of  temperature  only, 
whereas  Equation  (26)  is  a  function  of  time  and  temperature.     Based  on  this 
assumption,   the  stimulus  (I)   in  Equation  (22)  for  pain  (S)  can  be 
represented  by  the  damage  rate  (dl^/dt) ,     Thus  the  relationship  between  pain 
and  tissue  damage  can  be  expressed  according  to  the  Weber -Fechner  law 
(Equation  (22))  as: 

S  -  So  +  k  log  Q/Qo  (32) 

where  S    -  pain  intensity 
Sq  -  threshold  pain 
k    -  constant 
Q    -  dfl/dt  -  damage  rate 
Qq  -  damage  rate  at  pain  threshold 

As  a  first  attempt,  Hardy  et  al.  combined  Equations  (32)  and  (25)  to 

calculate  the  pain  intensity  and  its  corresponding  skin  surface  temperature 

using  the  rate  model  developed  by  Henriques.     The  calculated  values  were 

then  compared  with  the  experimental  values.     They  failed  at  pain 


intensities  above  17  JND  (8.5  dol,  see  Figure  12).     Later,  Stoll  and  Greene 
repeated  the  procedures  to  find  the  relationship  between  tissue  temperature 
at  80  /im,  damage  rate  and  pain  intensity.     This  time,  they  used  their  own 
damage  rate  data  (Figure  13)  instead  of  the  data  from  Henriques .  They 
successfully  obtained  the  correlation  shown  in  Figure  14.     It  should  be 
noted  that  the  JND  and  temperature  in  the  abscissa  of  Figure  14  are  not 
corresponding  to  each  other.     The  abscissa  actually  represents  two  axes 
overlapping  each  other.     That  is,  the  two  curves  in  the  figure  are  the 
projections  of  a  single  curve  in  space.     This  relationship  is  also 
tabulated  in  Table  13.     Since  it  has  been  shown  that  threshold  blistering 
occurs  within  24  hours  when      -  1,  and,  the  damage  rate  corresponding  to  a 
constant  tissue  temperature  or  pain  intensity  remains  constant,  then  the 
time  to  blistering  at  a  constant  tissue  temperature  or  pain  intensity  can 
be  predicted  by  the  reciprocal  of  the  damage  rate  at  which  the  temperature 
or  pain  intensity  is  maintained  (Table  13). 


Figure  14      Correlation  of  Pain  Intensity  with  Damage  Rate  and  Tissue 
Temperature  at  80  pm.     From  Stoll  and  Greene  [12] 


Table  13      Relationship  Between  Pain,  Tissue  Temperature  and  Damage  Rate. 
From  Stoll  and  Greene  [12] 


Observed  Tissue  Damage  Rate  Predicted  Time 

Pain  Intensity        Temperature  Q  H/Qq  to  Threshold 

at  80  fim  Blister  when  Q-1 

(JND)  ("C)  (sec-^)  (sec) 


1 

45 

0.00033 

1.0 

3300.0 

2 

0.00047 

1.4 

2130.0 

4 

46 

0.00097 

2.9 

1030.0 

8 

0.00420 

12.4 

228.0 

10 

48 

0.00840 

25.4 

118.0 

12 

0.01800 

54.5 

56.0 

14 

50 

0.03600 

109.0 

28.0 

16 

0.08300 

252.0 

12.0 

18 

53 

0.15000 

456.0 

6.7 

20 

0.30000 

912.0 

3.3 

21 

56 

0.47000 

1430.0 

2.1 

7.       COMPUTER  MODELS  TO  DETEEIMINE  THE  THERMAL  INJURY  FROM  EXPERIMENTAL  DATA 

Since  the  availability  of  adequate  information  on  the  thermal 
tolerances  and  damage  of  living  tissue  due  to  heating,  computer  models  have 
been  developed  for  the  analysis  of  experimental  data.     In  1960,  Stoll  [23] 
first  developed  an  algorithm  to  evaluate  the  measurements  of  the  changes  in 
surface  temperature  with  time  throughout  heating  and  cooling  of  the  skin  by 
thermal  irradiation.     From  these  measured  surface  temperatures  (Tg) ,  Stoll 
calculated  the  tissue  temperature  at  the  basal  layer  (approximately  80 
deep,  where  the  base  of  the  blister  forms)  using  Equation  (8)=  Typical 
measured  temperature  histories  at  the  skin  surface  and  those  computed  at  80 
/xra  are  shown  in  Figure  15  for  two  irradiances  with  durations  just 
sufficient  to  produce  minimal  blistering  within  24  hours.     The  damage  rate 
(dfl/dt)  was  then  determined  by  using  the  radiative  heating  curve  in  Figure 
13.     This  curve  is  represented  by: 

log  ^  -  slope  X  T^  +  log  (intercept)  (33) 

where        slope  -  0.42      for      T^  <  50'C 
-  0.16      for      Tx  >  50'C 
Tx  -  tissue  temperature  at  80  /im,  ®C. 

Figure  16  shows  the  plots  of  the  damage  rate  vs.  the  time  at  which  the 

basal  layer  temperature  occurred  for  the  two  irradiances.    To  calculate  the 

total  damage  (O)  including  the  damage  during  heating  as  well  as  cooling, 

Stoll  integrated  the  curves  shown  in  Figure  16  over  the  periods  the 

temperatures  at  the  basal  layer  remained  above  44*C  (the  injurious  level) . 

Applying  the  trapezoidal  rule  for  integration,  the  total  damage  can  be 

expressed  as: 
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1 


n-;  —  dt  +  J  — 


p  dn  r  dn 


dt 


i  dt  p  dt 


-  1^  [n   +  2Q.  +       +  . . .  +  2n 

2        o  1  2 


(N-1) 


(34) 


where    i    -  time  at  which  the  injurious  temperature  level  (44°C)  is 


reached,  sec 
p    -  time  at  which  the  exposure  is  removed,  sec 
ij«  -  time  at  which  temperature  returns  to  44'C,  sec 
At  -  time  interval  for  numerical  integration,  sec 
Qj  -  damage  rate  dfi/dt  at  the  j^^  grid  point,  sec"^ 
N    -  total  number  of  grid  points  including  the  heating  and  cooling 


The  above  procedure  developed  by  Stoll  [23]  requires  direct  measure- 
ment of  the  skin  surface  temperature  during  irradiation.    Weaver  and  Stoll 
[10]  modified  the  method  to  eliminate  the  measurements  of  skin  temperature. 
Instead,  they  measured  the  radiation  intensity  and  time  (a  much  easier  task 
as  mentioned  by  them)  and  calculated  the  basal  layer  temperature  using 
Equation  (11) .     They  also  considered  the  variations  of  thermal  conductivity 
of  the  skin  with  the  temperature  of  the  basal  layer.     By  assuming  pc  =  1 
cal/cm^-°C  so  that  a  »  k,  the  conductivity  solved  by  them  has  been 
presented  in  Figure  8  in  Chapter  4.     Since  little  variation  was  observed  at 
temperatures  above  44'C,  they  assumed  a  constant  value  for  the  thermal 
conductivity  during  the  heating  phases,  that  is  the  average  value  for 
temperatures  between  44 and  the  peak  temperature.    During  the  cooling 
phases,  it  was  found  that  a  constant  conductivity  failed  to  represent  the 
temperature  drop.    Weaver  and  Stoll,  then,  established  an  algorithm  to 
estimate  the  conductivity  during  cooling  at  each  temperature  point  corre- 
sponding to  each  of  the  grid  points  on  the  time  scale  after  the  exposure 
was  removed.    This  algorithm  is: 


phases 


y  k j  -  y  k^_^  -  10-^.     j  -  2,  3,   ...  n  (35) 

where    k]^  -  The  average  conductivity  during  cooling  for  temperatures 
greater  than  44''C  and  less  than  the  peak  temperature, 
cal/cm  sec*C 

n    -  number  of  grid  points  after  the  peak  temperature  point 
Thus,  the  square  root  of  the  conductivity  at  the  j^^  time- temperature  point 
was  estimated  by  subtracting  a  constant  value  of  10'^  (cal/cm  sec'C)^'''^ 
from  that  at  the  previous  point.     This  procedure  was  repeated  for  each 
successive  point  until  the  point  at  which  the  tissue  temperature  at  80  /xm 
deep  decreases  back  to  44°C,     Figure  17  shows  the  plots  of  the  theoretical 
conductivity  as  predicted  by  this  method.     Some  of  the  experimental  conduc- 
tivity curves  from  Figure  8  are  also  included  for  comparison.  This 
theoretical  function  of  conductivity  was  found  to  be  capable  of  repre- 
senting the  experimental  time- temperature  histories  at  various  levels  of 
radiation.     A  comparison  of  the  experimental  and  the  calculated  theoretical 
time -temperature  histories  for  a  level  of  radiation  0.4  cal/cm^  sec  is 
given  in  Figure  18. 

Once  the  time- temperature  history  at  the  basal  layer  was  obtained, 
Weaver  and  Stoll  computed  the  damage  rate  by  using  Equation  (25) .  The 
constant  P  and  the  activation  energy  divided  by  the  gas  constant  AE/R  were 
determined  by  using  Equation  (33)  or  Figure  13.    Values  of  dQ/dt  and  T^ 
from  Equation  (33)  were  substituted  into  Equation  (25)  which  yielded 
solutions  for  P  and  AE/R.     These  values  are: 


P    -  2.185  X  lO''-^^  sec""*- 


AE/R  -  93534.9  K 
and 

P  -  1.8230  X  10^^  sec'"*- 
AE/R  -  39109.8  K 


for    44'C  <  Tx  <  50'C 


for  Tx  >  50'C 
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Figure  17    Experimental  and  Theoretical 

Thermal  Conductivity  of  Skin 

versus  Tissue  Temperature.  ^ 
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Figure  18    Experimental  and  Theoretical 
Time -Temperature  History  at 
the  Basal  Layer  for  an 
Irradiance  of  0.4  cal/cm^sec. 
From  Weaver  and  Stoll  [10].  ssf 
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Similarly,  upon  integration  of  the  damage  rate  by  the  trapezoidal  rule 
gives  the  total  damage. 

The  exposure  times  (r)  required  to  produce  minimal  blistering  within 
24  hours  for  irradiation  between  0.1  and  O.A  cal/cm^sec  have  been  measured 
[12]  and  the  changes  in  thermal  conductivity  with  temperature  during  the 
heating  and  cooling  phases  for  this  range  of  irradiation  have  also  been 
predicted  [10].     In  order  to  obtain  the  same  information  for  irradiation 
above  0.4  cal/cm^sec,  Weaver  and  Stoll  plotted  the  exposure  time  (r) ,  the 
thermal  conductivity  during  heating,   and  that  during  cooling  against  the 
radiation  with  extrapolations  to  the  higher  levels  (Figures  19  and  20) . 
The  effective  radiation  indicated  in  Figures  19  and  20  is  the  absorbed 
radiation  by  the  skin.     In  Figure  20  the  value  of  the  thermal  conductivity 
during  heating  for  any  given  radiation  level  is  the  constant  value  used  to 
compute  the  temperature  rise  at  every  time- temperature  point.     However,  the 
conductivity  during  cooling  is  k]^,  the  conductivity  at  the  first  time- 
temperature  point  after  the  peak  temperature.     From  the  information  given 
in  Figures  19  and  20,  the  tissue  temperature  at  80  /im  deep  when  t  -  r  (peak 
T^) ,  the  total  time  from  t  -  0  to  the  time  at  which  the  tissue  temperature 
falls  below  44*C,  the  damage  during  heating  (fin)*  damage  during  cooling 

(fi^) .  the  total  damage  (Q-p  -       +  Oc)  ^  radiation  level  of  3.0 

cal/cm^sec  can  be  calculated  by  means  of  the  numerical  procedure  developed 
by  Weaver  and  Stoll.     These  results  are  summarized  in  Table  14  in  which  Q 
is  the  effective  radiation;  At  is  the  time  interval  used  to  calculate  0; 
and  ki  and  A2  are  the  square  root  of  the  conductivity  shown  in  Figure  20 
during  heating  and  cooling  respectively. 
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Figure  19      Reciprocal  of  r  versus  Radiation  Level.     From  Weaver  and  Stoll 
[10] 
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Figure  20      Thermal  Conductivity  versus  Radiation  Level.     From  Weaver  and 
Stoll  [10] 


As  seen,  n-p  =  1  for  all  the  exposures  shown  in  Table  14.  Therefore, 
they  are  the  threshold  exposures  just  sufficient  to  cause  minimal 
blistering.     Also,  it  is  noted  that  a  larger  portion  of  the  tissue  damage 
is  done  after  the  exposure  ceased  at  the  higher  radiation  levels.     In  this 
case,   the  total  amount  of  heat  required  to  produce  minimal  blistering  is 
lesser  at  the  higher  radiation  levels  than  at  the  lower  levels.     The  damage 
incurred  during  cooling  can  be  minimized  by  prompt  application  of  ice  to 
the  exposed  area. 


Table  14        Characteristics  of  Time -Temperature  Histories  of  Skin  Exposed 
to  Thermal  Radiation  Just  Sufficient  to  Cause  Minimal 
Blistering.     From  Weaver  and  Stoll  [10]. 
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In  the  1970' s,  skin  damage  due  to  laser  irradiation  became  a  topical 
subject.     Consequently,  computer  models  that  are  suitable  for  analyzing  the 
complex  conditions  of  burn  injury  produced  by  a  laser  were  developed.  These 
models  are  capable  of  simulating  the  effects  of  blood  flow,  metabolism,  and 
temporal  and  multi- dimensional  spatial  variations  in  heating  [18,  24].  In 
the  present  review,  laser  irradiation  of  tissue  will  not  be  covered. 


8.  SUMMARY 

The  previous  sections  have  described  the  mathematical  and  physical 
analysis  of  the  physiological  response  of  skin  to  burns.     Methods  to 
predict  thermal  damage  to  tissue  have  also  been  discussed. 

Under  normal  conditions,   the  skin  temperature  is  about  33°C.     When  the 
skin  is  exposed  to  a  heat  flux,  the  skin  temperature  increases.  The 
increase  of  skin  temperature  can  be  described  by  the  standard  one- 
dimensional  heat  conduction  equation.     When  the  temperature  of  the  skin  is 
elevated  to  about  44*'C,   the  pain  receptors  in  the  skin  will  give  warning  of 
the  potentially  harmful  conditions.     At  about  the  same  temperature,  tissue 
damage  begins.     The  rate  of  damage  can  be  represented  by  a  rate  process 
equation  which  is  a  function  of  the  skin  temperature  only.     An  arbitrary 
function  Q  is  used  in  this  equation  to  describe  the  amount  of  damage.  When 
Q  -  1,  the  threshold  of  second  degree  burn  occurs.     Second  degree  burn  is 
defined  as  the  damage  that  Indicates  complete  destruction  of  the  epidermis 
without  a  significant  degree  of  injury  to  the  dermis.    At  this  threshold 
burn,  the  attachment  between  epidermis  and  dermis  is  impaired.     The  amount 
of  edema  fluid  that  accumulated  under  the  epidermis  is  usually  sufficient 
to  form  a  blister.     Thus  blistering  normally  indicates  destruction  of 
epidermis.    When  Q  >  1,  damages  are  irreversible. 

Like  the  rate  of  tissue  damage,  pain  depends  only  on  the  temperature 
to  which  the  skin  is  raised.     Maximum  perceptible  pain  occurs  at  around  60 
to  yO^C.    While  the  amount  of  damage  (degree  bums)  and  blistering  are 
dependent  on  both  the  temperature  and  time.     The  temperature  and  time 
required  to  produce  a  certian  degree  of  damage  are  inversely  related  and 
can  be  described  by  a  simple  exponential  function.     In  general,  when  the 


temperature  of  the  skin  is  maintained  at  44 (the  minimum  injurious 
level) ,  an  exposure  of  about  6  hours  is  required  before  the  occurrence  of  a 
second  degree  burn.     At  temperatures  above  70"C,  a  second  degree  burn  can 
be  produced  in  less  than  1  second. 

The  occurrence  of  pain  and  tissue  damage  during  a  burn  are  mainly  due 
to  the  thermal  denaturation  of  cellular  proteins  which  have  an  activation 
energy  of  roughly  150,000  cal/mole. 
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